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Abstract 

We study classical and quantum dynamics on the Euclidean Taub-NUT geometry coupled 
to an abelian gauge field with self-dual curvature and show that, even though Taub-NUT 
has neither bounded orbits nor quantum bound states, the magnetic binding via the gauge 
field produces both. The conserved Runge-Lenz vector of Taub-NUT dynamics survives, in 
a modified form, in the gauged model and allows for an essentially algebraic computation of 
classical trajectories and energies of quantum bound states. We also compute scattering 
cross sections and find a surprising electric-magnetic duality. Finally, we exhibit the 
dynamical symmetry behind the conserved Runge-Lenz and angular momentum vectors 
in terms of a twistorial formulation of phase space. 


1 Introduction 

The Euclidean Taub-NUT (TN) geometry has been studied extensively and from several dif¬ 
ferent points of view. It is interesting as a particularly simple example of a gravitational 
instanton, it can be viewed as a Kaluza-Klein geometrisation of the Dirac monopole and it 
arises in the context of monopole moduli spaces, either directly, or, in its ‘negative mass’ form, 
as an asymptotic limit. In each of these contexts it plays a role akin to that of the hydrogen 
atom in atomic physics, both in the general sense of being the simplest example and in the 
technical sense of sharing simplifying features with the hydrogen atom, like a Runge-Lenz type 
conserved quantity. 

The four-dimensional Maxwell equations on TN space have a simple source-free solution which 
is intimately connected to the TN geometry. This was hrst pointed out by Pope niq who went 
on to show that that the index of the Dirac operator minimally coupled to this Maxwell held 
is non-trivial. The index and the properties of zero-modes of this gauged Dirac operator were 
recently studied in detail in our paper [3] from which the current paper evolved. 

The Maxwell held hrst considered by Pope has played an important role in various contexts. Its 
held strength is a harmonic two-form which is self-dual for an appropriate choice of orientation 
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and also square-integrable. It is exact, with a globally defined gauge potential which is, however, 
not square-integrable. In other words, the harmonic two-form generates the non-trivial L^- 
cohomology in the middle dimension of TN space. This is the reason why it was important in 
tests of S-duality on monopole moduli spaces HE]. 

One can relate the self-dual two-form directly to the TN geometry by noting that, with a 
suitable normalisation, it is the Poincare dual of the CP^ which compactifies TN to CP^ [B]. 
More generally, one can understand the L^-cohomology of TN and its multi-centre generalisation 
in terms of the ordinary cohomology of a suitable compactification [7j. 

For all of these reasons, it is not surprising that the inclusion of the self-dual gauge field in 
the dynamics on TN space turns out to be mathematically natural. In this paper, we consider 
the geodesic motion on TN coupled to the self-dual gauge field as our model for the classical 
dynamics and the Dirac and Laplace operators on TN, also minimally coupled to the gauge 
field, as quantum models. We show that all the interesting algebraic features of ordinary TN 
dynamics carry over to the gauged case, and that, moreover, bounded motions and quantum 
bound states, neither of which are possible on (regular, ‘positive mass’) TN alone, occur in the 
gauged dynamics. 

In general qualitative and semi-quantitative terms, the classical and quantum dynamics on TN 
space is strikingly similar to the Kepler problem and the non-relativistic hydrogen atom. The 
basic reason why the four-dimensional TN geometry can model the three-dimensional motion of 
a charged particle is the fact that TN is a Kaluza-Klein geometrisation of the Dirac monopole in 
three dimensions [HE]. The more detailed similarities are related to the conserved Runge-Lenz 
vector in both cases. In this sense, one can view our inclusion of the magnetic field in the 
TN dynamics as analogous to the inclusion of a magnetic monopole field in the much studied 
extension of the Kepler problem to the MICZ Kepler problem mm- 

From the Kaluza-Klein point of view, the magnetic field of the monopole is already encoded 
in the geometry, and the name ‘magnetic field’ for an additional abelian gauge field on TN is 
potentially confusing. We adopt it here because it is justified from the four-dimensional point 
of view. As we shall see, the magnetic field leads to magnetic binding akin to that responsible 
for Landau levels in planar systems. In fact, in a limit where TN becomes fiat Euclidean 
four-space, the four-dimensional magnetic field is constant, and the bound states that we find 
become ordinary Landau levels. This picture of magnetic binding also provides a qualitative 
explanation of the index found by Pope and of the form of the zero-modes discussed in j2]. 

We have organised our presentation to proceed from the most direct to more abstract treatments 
of gauged TN dynamics. We begin in Sect. 2 with a brief general discussion of how a magnetic 
field on a two-dimensional Riemannian manifold can produce bound states. In Sect. 3, we collect 
conventions for describing the TN geometry and the associated Dirac and Laplace operators 
coupled to a magnetic field. We turn to the classical dynamics in Sect. 4, discuss the conserved 
Runge-Lenz and angular momentum vectors of the gauged geodesic motion on TN and describe 


2 


the classical trajectories. Sect. 5 contains a direct solntion of the eigenvalue problem for the 
gauged Laplace operator on TN space through separation of variables. We exhibit the promised 
bound states, give their energies and degeneracies and compute scattering cross sections. In 
Sect. 6 , we solve the quantum problem algebraically, using a quantum version of the Runge-Lenz 
vector. In Sect. 7, we exhibit the symmetry underlying the conservation of angular momentum 
and the Runge-Lenz vectors from a twistorial description of phase space. Our hnal Sect. 8 
contains a brief discussion, our conclusions and an outlook onto open problems. 

2 A toy model: motion on a surface with magnetic field 

We can gain a qualitative understanding of bound states on TN coupled to a Maxwell held by 
considering a two-dimensional model, consisting of a two-dimensional manifold with metric and 
magnetic held. We will encounter a manifold and metric of the same kind in our study of TN 
as a geodesic submanifold, and the magnetic held as the restriction of the Maxwell held to the 
geodesic submanifold. However, here we study the two-dimensional model in its own right. 

Consider a two-dimensional manifold diheomorphic to an open disk D with f/(l)-invariant 
metric of the form 

ds'^ = dR^ + c\R)dj\ ( 2 . 1 ) 

For consistency with our later discussion of the TN geometry we take the angular coordinate 7 
in the interval [ 0 , 47 r), so that dvrc is the length of a U{1) orbit. The radial coordinate R is the 
proper radial distance from the origin and has range [0, cxd), and we assume a form of c near 
R = 0 to ensure that the metric is smooth there. We are interested in two kinds of behaviour 
of the function c. 

The hrst case captures what happens in the regular TN geometry. The function c has the 
hnite range [0, L) for some positive real number L so that the length of the U{1) orbits remain 
bounded. Moreover we assume that c(0) = 0 and that c is strictly monotonic, so that one can 
picture the metric as being induced on a cigar-shaped surface of revolution in three-dimensional 
Euclidean space, as shown in Fig. [T} The qualitative behaviour of geodesics on such a surface is 
well know and follows from Clairaut’s relation. Generic geodesics spiral on the cigar. Geodesics 
spiralling towards the tip will be reflected at some point and spiral out. All geodesics ultimately 
move arbitrarily far away from the tip and there are no geodesics which remain in a region 
bounded by a hnite value of R. 

The second case captures what happens in the singular or ‘negative mass’ TN. The function 
c diverges at R = 0 , has the range (L, 00 ) and is monotonically decreasing. As an embedded 
surface, this is a funnel, with the opening at R = 0 and the tip at R = 00 as shown in 
Fig.B Generic geodesics again spiral on this surface, but now there are two kinds of behaviour. 
Geodesics which travel straight down the funnel or spiral only slowly may escape to R = cxd. 
However, geodesics travelling into the funnel with sufficiently high angular momentum relative 
to their speed will bounce back and remain inside a region bounded by some hnite value of R. 
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Figure 1: The cigar-shaped surface for positive L (left) and the funnel-shaped surface for 
negative L (right). 


We now return to the hrst case with monotonically increasing c G [0, L) and consider the 
inclusion of a magnetic held of a specihc type given by the two-form 


B = d 


pc , \ P , , 


( 2 . 2 ) 


for some real constant p which controls the strength of the magnetic held, and is proportional 
to its hux: 


— / B = p. 
271 ' ^ 


(2.3) 


' D 


The Lagrangian governing the motion of particle on the surface with metric (2.1), minimally 
coupled to the gauge potential for B is, for a suitably chosen mass parameter. 


£ = -(i?2 + cV 


pc 


7- 


With the momenta conjugate to R and 7 

dC 1 . 
^ dR 2 


the Hamiltonian is 


dC In pc^ 


q pc 


«-PH+(-+2i2 


(2.4) 


(2.5) 


(2.6) 


Since q is conserved and p constant, this is ehectively the Hamiltonian for one-dimensional 
motion on the half-line in the potential 




(2.7) 


We would like to know if there are bounded trajectories in the potential ( |2.7 ). As a potential, 
W should be viewed as a function of i?, but with our assumption that dc/dR > 0 we can study 
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Figure 2: Plots of the potential (2.7) for L 
q = 3 and p = 8 (right). 



1 and c G [0,1) for g = 3 and p = —8 (left) and 


its minima by looking at PF as a function of c. It is easy to check that W{c) has a unique 
minimum at Cm > 0 satisfying 

^2 On 

( 2 . 8 ) 


L2 


2q 

p 


However, for Cm to be in the range [0, L) we require 

|9| < 


(2.9) 


and this is a necessary and sufficient condition for W to have a minimum. The value at the 
minimum is 

0 if pg < 0 


W{Cm) = 


^ ifp^>0. 


( 2 . 10 ) 


The qualitative form of the potential is similar in the two cases. We also note that the inequality 

( 2 . 11 ) 


(2.10) implies the bound 


L'^W > 


0 if pg < 0 
2pg if pg > 0, 

which will play an important role in our discussion. 

We conclude that a magnetic field on a cigar-shaped surface can lead to bounded trajectories 
even though the geometry of the cigar does not support any bounded geodesics. In the presence 
of the magnetic field (2.2), all trajectories with angular momentum g of magnitude less than 
|p/2| remain in a bounded region. 
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3 Classical and quantum geometry of the Taub-NUT space 

3.1 The Taub-NUT geometry 


We use the conventions of [3] for the description of the TN geometry, and refer the reader to 
that paper for details. Here we summarise essential background and notation. The TN space 
Mtn is a non-compact four-dimensional Riemannian manifold with a self-dual Riemann tensor 
and U{2) isometry. As a manifold it is diffeomorphic to and the metric can be written as 

+ c^crl, (3.1) 


where aj are the left-invariant one-forms on SU{2) ~ S^. They are dehned in terms of h G 
SU{2) via 

h~^dh = aiti (72^2 + (3.2) 

where we use su{2) generators which are given in terms of the Pauli matrices tj as tj = 

The quantities a,b,c,f are functions of a radial coordinate r transverse to the SU{2) orbits. 
The self-duality of the Riemann curvature is then equivalent to 

^^ = (&-c)^ - a^+cycl. , (3.3) 

and the solution which gives rise to the TN metric is 


= b = c = / 




(3.4) 


Here e and L are parameters which are required to be positive for a smooth metric. As discussed 
in na, the parameter e is relevant for taking the Landau limit of our model, and can be used 
to introduce time-dependence in an interesting way. However, we shall set e = 1 for most of 
the text and work with the potential 

U = l + -, (3.5) 

r 

unless stated otherwise. The TN space with L < 0 arises as the asymptotic form of the two 
monopole moduli space na but has a singularity at r = —L. We will assume L > 0 in the 
following. 

For the discussion of symmetries and the definition of natural differential operators associated 
to the TN geometry, we require the vector fields Xj which are dual to the left-invariant forms 
CTj (so = 6ij) and which generate right-actions on h G SU{2), 

Xj : h ^ htj. (3.6) 

We will also need the vector fields Zj which generate the left-action, 

Zj h —y — tjh. (3-7) 
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It follows from these definitions that, as vector helds on S^, 

[Xi^Xj] = eijkXk, [Zi, Zj] = eijkZk, [Xj, = 0, (3.8) 

where summation over repeated indices is assumed. By dehnition, the Lie derivative of the 
one-forms cTj with respect to the Zj vanishes. 

There are several natural ways to parametrise SU{2) and hence to write the one-forms and 
vector helds on it. For the purposes of this paper we need both a parametrisation in terms of 
Euler angles and one in terms of two complex numbers. Continuing with the conventions of 
[3] we therefore introduce complex numbers Zi,Z 2 , constrained to satisfy -f |z 2 p = 1, and 
write elements of SU{2) as 



V 2:2 Zi J 

(3.9) 

In some calculation it is more convenient to work with Euler angles /3 G [0, tt], 

a G [0, 27 r ),7 G 

[ 0 , dvr) dehned via 

— 4 (ck+'y) ^ —T') ’ ^ 

zi = e (.Qg _^ ^ g 2 sin — . 

(3.10) 

Then the left-invariant one-forms can be expressed in terms of Euler angles as 



(Ti = sin 7(7/3 — COS 7 sin 

(72 = cos 7 ( 7 / 3 -|-sin 7 sin/3(iQ;, 

(73 = ( 77 -f cos/3(7a;, 

(3.11) 

and the dual vector helds 

are 



Xi = cot /3 cos + sin - 

smp 

X 2 = — cot (5 sin-f cos pdn H- %,da, 

sm /3 



^ 3 = d^. 

(3.12) 

In terms of complex coordinates we have 



Xi — -( 2 : 2^1 ~ ~ ^ 2^1 + ^ 182 ), 

X 2 = ^(-72^1 -F Zid 2 - Z 2 di -7 Zid 2 ), 

X 3 = -(—2;ic7i — ^ 2^2 + Zidi + Z 2 B 2 ), 

(3.13) 

and 

Zi = ^{-Z 2 di - zid 2 Z 2 di + 01 ^ 2 ), 

Z 2 = -(—^ 2^1 + ^ 1^2 ~ ^281 -f ^ 1 ^ 2 ), 

% — — 

Z 3 = — Z 282 — Zi 8 i + ^ 282 ). 

(3.14) 
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It follows from general arguments, but can also be checked explicitly, that the Laplace operator 
on ~ SU{2) can be written as 


A53 = 


(3.16) 


The magnetic held on TN referred to in the title and discussed in the Introduction can be 
written as the exterior derivative of the abelian gauge held 

2 

T) r 

(3.16) 


pc p r 

= -—r^cTa, 


2L2 " 2r + L' 

where p is a real parameter which will play an important role in our discussions. The magnetic 
held 


jr = dA=f^ 


cdc A cTa + —ct 2 A ai 


(3.17) 


is (up to a multiplicative constant) the unique harmonic, normalisable and t/( 2 )-invariant two- 
form on TN. 

Note that the gauge we have chosen preserves the U{2) symmetry of the TN geometry. However, 
the magnetic held inevitably breaks the discrete symmetry 


which maps 


/3i—^-TT — /3, a ^ 71 + a, 71 —)■ — 7 , 


(Ji H->• (Ti, (72 —(J 2 , (7a —(7a, 


(3.18) 


(3.19) 


and therefore preserves the metric (3.1) but neither the gauge held (3.16) nor its curvature. As 


we will see, this has interesting consequences in the dynamics. 

TN without the point r = 0 (the NUT) is a circle bundle over \ {0}. As reviewed in 
for each direction {(3, a) in the base there is a geodesic submanifold parametrised by ( 7 , 7 ). 
Each of these geodesic submanifolds is of the general cigar-shape of our toy model in Sect. 2, 


and the hux of the magnetic held (3.17) through this submanifold is 27rp. Thinking of TN as 


a two-sphere’s worth of such cigar-shaped surfaces threaded by magnetic hux will prove very 
helpful for a qualitative understanding of our results. 


3.2 The gauged Dirac and Laplace operators 


The most fundamental operator associated to the metric (3.1) and the connection (3.16) is the 


Dirac operator on TN minimally coupled to the connection. As shown in HD for the ungauged 
case, the spectrum of the Dirac operator and the Laplace operator are closely related. The 
arguments in [T3] are essentially a rehection of an underlying supersymmetry. It is not difficult 
to adapt them to the gauged case, as we shall now show. 

The gauged Dirac operator on TN has the following form [3] 


lPr,= 


0 


p 

0 


(3.20) 









where 


Lf- 


= 

" V 

w(' 


-dr 

-dr 




1 V 

r^lL^ 2r2V 


+ ^(nXi + T2X2) 


(3.21) 


As shown in DEI and elaborated in j3], the kernel of is trivial but that of Tp has dimension 
|[|p|]([|p|] + 1)5 where [x] is the largest integer strictly smaller than the positive real number x. 
It follows that 

H. = TpT^ (3.22) 

is a strictly positive operator, and that 

H+ = T^T, (3.23) 

is positive, but not strictly positive. Therefore, we can dehne the unitary operator 

u = (3.24) 

with inverse U~^ = T^/^/HZ, and use it to relate and via 

= U-^H_U. (3.25) 


It follows from this unitary equivalence that, as in the ungauged case [H], has the same 
spectrum as apart from the zero eigenvalue of . In other words if T is a two-component 
eigenspinor of H_ with eigenvalue E then is an eigenspinor of with the same eigen¬ 

value. 


Combining these results, we obtain eigenstates of the Dirac operator with non-zero eigenvalues 
from eigenstates of H_ as follows. Suppose that T is an eigenstate of H_ with eigenvalue > 0. 
Then 



(3.26) 


Inserting the expressions for the TN prohle functions, we hud 


H = T =_ — 

^P-^P p2V- 


1/ 


dr{r‘^dr) - [(Xi ti)^ {X 2 -h ^ 2 )^] ~ 7^ ("^3 + 7^ + ^3 ) • (3.27) 


L2 


ip 


2V 


It is convenient to change gauge by observing from (3.6) that 

Xih~^ = —tih~^ 


(3.28) 


and that therefore, as an operator identity, 

h{Xi + ti)h-^ = h{Xih-^) + hh-^Xi + hUh-^ = Xi. (3.29) 
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Employing this, we obtain 


(3.30) 


hH.h-^ = HpTo, 

where tq is the 2 x 2 identity matrix and 


Hp = 


1 


r‘^V 


dr{r‘^dr) 




(3.31) 


is the Laplace operator associated to the metric (3.1) and minimally coupled to the gauge held 


(3.16). This is the operator whose spectrum we shall study in the remainder of this paper. 


4 Dynamical symmetries in classical Taub-NUT dynamics 

4.1 Canonical procedure 

We now turn our attention to the classical dynamics in TN in the gauged case. We discuss the 
conserved angular momentum and Runge-Lenz vectors, and use them to describe the classical 
trajectories. Our treatment is an extension of the discussion in [I3] and [I5] of the (ungauged) 
motion on TN space. 


As reviewed in Sect. 2, the TN space (3.1) can be parametrised by a coordinate r and the Euler 


angles a,/3 and 7 . In these coordinates the Lagrangian for geodesic motion on TN takes the 
form 

+ + ( 4 . 1 ) 

where Ui are the components of the body hxed angular velocity. 


oji = siny/3 — cosy siu/S d, 
0 J 2 = cosy/3 + sin y sin/3 d, 
1^3 = y + cos/9d. 


(4.2) 


and we have chosen an overall factor of 1/4 for convenience. 

Inserting the angular velocities in the Lagrangian and recalling that a = b = r\fV^ c = 
L/yr,/ = -Vr we obtain 




R(f 2 + ^2 g .^2 ^^ 2 ) ^ LV-^(y + cos Paf 

In terms of cartesian coordinates 


r = {xi,X 2 -iXz) = (r sin/3 cos a, r sin/3sina, r cos/9). 


(4.3) 


(4.4) 


the Lagrangian takes the more familiar form 


C = + LV-i(y + A ■ T)2), 


(4.5) 
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where A is a gauge potential for the Dirac monopole 

X3X2 , X3X1 


Ai — 


r{r‘^ — x|) 


, ^2 — 


which satishes 


as well as 


r(r 2 - xl) 


Xk 


, ^3 — 0 , 


dlAra - draAi = Sklm^ for r 7^ 0, 


A - dr = cos (3 da. 


(4.6) 


(4.7) 


(4.8) 


We now minimally couple the motion on TN to the gauge potential (3.16) via the Lagrangian 

+ ^(7 + ^(7 + A -^. 

Clearly, the momentum q conjugate to the cyclic coordinate 7 , 

dCp • p 

is conserved. The canonical momentum vf conjugate to r’is 

A 

2V 

where 


47(7 + ^-^^“ ^A = p + gA, 


dr 


p = -Vr 
2 

is called the mechanical momentum |13j . 

The canonical symplectic structures on the phase space T*Mxn, 

dxi A dni + d'y A dq, 


(4.9) 

(4.10) 

(4.11) 

(4.12) 


(4.13) 


is invariant under the 17(1) action which maps 7 —)■ 7 + 5. The moment map for this action is 
the charge q, viewed as map —>■ M, and the symplectic quotient by this 17(1) action 


= T*M™//t/(l) 


(4.14) 


is, by dehnition, the pre-image of any real constant under the map q divided by the 17(1) action. 
The position vector r and the canonical momentum vector vf provide natural coordinates in 
terms of which the symplectic structure on A4q takes the form 


q 

u = dxi A d-Ki = dxi A dpi -7 ^eunXndxi A dxi. 

Ir^ 


The associated Poisson brackets are 


dAdB dA dB 


(4.15) 


(4.16) 
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(4.17) 


SO that the mechanical momentum p = n — qA satishes 


Xk 




where / is any function of r. 

We now rewrite the Lagrangian in terms of p and q, 


p 


and perform the Legendre transformation to obtain the gauged Hamiltonian, 

Hp = f-n + 'jq-Cp 

1 n 2 Tr pq p^ 

L2 L2 ALW 
= H + AH. 

Here H is the Hamiltonian for p = 0 and AH the contribution of the gauge potential: 

tr _ ^ m2 , ATT _pq , P^ 

y\^ +1^2^^ ALW 


(4.18) 


(4.19) 


(4.20) 


Recalling that the prohle function c appearing in the TN metric is c = L/\/V, we note that 

1 _ /« r4rN2 




> |? + ^V 

- \c 2H) 


> 


0 if pg < 0 
2pq if pq > 0. 


(4.21) 


For the last step we observed that the second term in the hrst line is the potential (2.7) of the 
toy model of Sect. 2, and used the bound (2.11). 


There is a conserved angular momentum ra of H given by 

J = r X p + gr, 


which, by virtue of (4.17), satishes the relations 

{JkiPl} ^klmPm- 

It also follows that 

{Jki Jl\ ^klmJm- 


(4.22) 


(4.23) 


(4.24) 


Relation (4.23) can be employed to check that J Poisson commutes with the Hamiltonian H. 
Since AH is spherically symmetric, J also commutes with Hp = H + AH. 
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In their study of the geodesic motion on the negative mass TN space in [ 13 ], Gibbons and 
Manton showed that there is a conserved vector quantity analogous to the Runge-Lenz vector 
of the Kepler problem which takes the form 


M = px J- ^ [L'^H - 2q^) . 
ZJj 


(4.25) 


One checks that it satishes 


\_Jki 


(4.26) 


and commutes with the TN Hamiltonian H for any value (positive or negative) of L. However, 
it fails to commute with our gauged Hamiltonian Hp since 


{AH, M,} = - 

By trial and error we hud that the vector-valued function 


/ = 


9 —* 

p r 
8 LrV' 


(4.27) 


(4.28) 


satishes {Hp, f^} = [AH, M^}. Hence the components of the gauged Runge-Lenz vector 

MP = M-f = px J-^{L^Hp-2q^-pq), (4.29) 

commute with Hp. The Poisson brackets between the components of J and turn out to be 

r 1 / n\2 1 




p(9+f) -H, 


(4.30) 


We will study their Lie-algebraic interpretation in detail in Sects. 6 and 7. 


4.2 Classical trajectories 

The conserved quantities discussed above can be used to determine the classical trajectories 


on TN in the gauged situation, i.e., the solutions of the Euler-Lagrange equations of (4.9) or 


Hamilton equation of (4.19) with Poisson brackets (4.16). Considering hrst the simpler case 


where g = 0, we deduce from (4.22) and (4.29) that 


J-r = 0, = MP -7^= - -LEr, 


(4.31) 


where E denotes the (constant) value of Hp and J = |J|. The hrst and second equations show 
that the movement is in a plane orthogonal to J and that is in this plane. Using polar 
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coordinates (r, 0 ) to parametrise the plane with in the direction determined by 0 = 0 , we 
dednce, from the third equation, 


r = -. 

\Mp\ COS0 + ^LE 

This is the equation of a conic section. Finally taking into account the relation 


(4.32) 


n>l = - 1^) + \l^E\ (4.33) 

we obtain the following types of orbit: the conic section is an ellipse for L^E < a parabola 
for L'^E = Y and a hyperbola for L'^E > 


J 



Figure 3: The conic sections determined by the conserved vectors J and N. 


In the general case O' 7 ^ 0, the expression (4.22) implies 


J ■ r = q, 


(4.34) 


which shows that f lies on a cone whose axis of symmetry is along J and whose vertex is at the 
origin. The opening angle 29 G (0, vr) of the cone relative to direction of J as shown in Fig. 
is determined by 


cos 6 


M 

J ■ 


(4.35) 
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For q > 0, the cone is in the ‘positive’ half-space determined by J ■ F > 0, while g < 0 it is in 
the ‘negative’ half-space determined by J ■ F< 0 . 

Furthermore, the equations (4.22[) and (|4.29) imply 


j.MP = -^{L^E-2q^-pq), ■ f = - q^ - ^ {L'^E - 2q^ - pq) . (4.36) 

ZIj ZJj 


To interpret them, we dehne the vector 


N = qMP + — {L^E - 2q^ - pq) J. 

ZIj 


(4.37) 


As a linear combination of conserved vectors with conserved coefficients, this vector is also 

(4.38) 


conserved. In terms of this vector, the second equation in (4.36) is equivalent to 

N ■ f= q{j‘^ — g^). 


which shows that the motion is also in a plane perpendicular to the vector N. With the notation 
Z = |F X p| for the magnitude of the orbital angular momentum, we note 


so that 


= /2 + 


N ■ f = qV 


(4.39) 


(4.40) 


The classical trajectories in the case g 7 ^ 0 are thus intersections of the cone dehned by (4.34) 


and the plane dehned by (4.40). From classical geometry we know that these are ellipses 


(including the degenerate case of a point), parabolae or hyperbolae (including the degenerate 
case of a line). The nature of the orbit depends on the energy E and on the relative size of g 
and p; as we shall see, the details are quite subtle, combining the results from the toy model 
in Sect. 2 with lessons from the role of conic sections as trajectories in the standard Kepler 
problem. 

Focusing on the non-degenerate case / 7 ^ 0, we note that the sign of g determines both the 


direction of the cone (4.34) and the position of the plane (4.40) relative to the origin. If g > 0 


then the situation is as shown in Fig. with the cone in the positive half-space determined by 


J-F> 0 and the plane (4.40) displaced from the origin in the direction of A^. If g < 0 the cone 
is in the opposite half-space and the plane is displaced from the origin in the direction of —N. 
The nature of the intersection between them, however, is independent of the sign of g, and only 
depends on the angle between J and iV, see again Fig. 

A lengthy calculation shows that the squared norm of N is 


fE 


|7V|2 = [L^E - 2pq) , 


(4.41) 
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which is positive for all allowed values of the energy by virtue of (4.21). Since, from the hrst 

P 


equation in (4.36), 


N-J=—{L^E-2q^-pq), 


we deduce that 6 is determined by 


cos (5 = 


I L^E-2q‘^-pq 
L^E{L^E -2pq) 


(4.42) 


(4.43) 


In order to classify the orbits we also note that, from (4.35) and (4.39), sin^ = j or 


TT 


COS ( — — 6^) = —. 


(4.44) 


Elementary geometrical considerations in Fig. |^now show that 



1 

V 

N._ 


" ellipse 

and < 

6 = ^-9 

> -v^ orbit is < 

parabola 


^-9 <S <^ + 9 


hyperbola 


.s>-,+e 


^empty set. 


(4.45) 


We analyse each of those conditions in turn. Since the cosine function is strictly decreasing on 


the interval [0,7r], applying it to the inequalities in (4.45) reverses them. It will also be useful 


to observe that the energy bound (4.21) implies 


^ ^ LPE > 2q^ + pq. 


(4.46) 


For the short proof, one needs to distinguish the cases pg > 0 and pg < 0 and use |g| < |p/2|. 

For elliptic orbits, we require cos5 > cos(| — 9). Inserting the above relations, this condition 
gives 

/ ^ 0 and L^E - 2q^ - pq > L^/E{L‘^E - 2pq). (4.47) 

Since the right hand side is positive (assuming L > 0), we deduce that, for elliptic orbits. 


L^E > 2q^ + pq. 


(4.48) 


On the other hand, squaring both sides of (4.47), we deduce 

L^E<U + P" 


2J 


However, the inequalities (4.48) and (4.49) can only both be satished if 

2 . 

^ < 4’ 


(4.49) 


(4.50) 
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which is precisely the condition (2.9) derived in the toy model in Sect. 2 . Since, by (4.46), the 


condition (4.50) is sufficient for (4.48) to hold, we deduce that elliptical orbits occur iff p 7 ^ 0, 


the charge q satisfy (4.50) and the energy satishes ( |4.49| ) As an aside we note that elliptical 
orbits are possible in the case L < 0 even when p = 0 (as discussed in [I3]). 


Returning to general p and positive L, the analysis of the conditions (4.45) for the parabolic 


and hyperbolic cases along the lines of the discussion of elliptical orbits is now straightforward. 
We skip most details, but point out that, in the hyperbolic case, the trigonometric identity 


cos (1 + 9) = — cos (f — 6 *) applied to (|4.45[) implies the condition 


/fyO and \L'^E — 2q^ — pq\ < L\/E{L‘^E — 2pq), 


(4.51) 


which (for positive L) is equivalent to 


L‘‘E> (, + | 


(4.62) 


but does not require any restrictions on p and q. 

We summarise the dependence of the orbits on the energy E and the charge q as follows: 


1 

\p^0, L^E<(q+l)^ 1 

1 1 

I ellipse 1 

1 fy 0 and < 

Wo, =(?+!)" 

> orbit is < 

parabola / 

1 

>(, + !)" J 

1 1 

[ hyperbola. J 


(4.53) 


5 Gauged Taub-NUT quantum mechanics 

5.1 Canonical qnantisation 

In this paper we set h = 1 when discussing quantum mechanics. With this convention, the 
canonical quantisation procedure of amounts to replacing 


^ .0 

TT q -^o^, 

or 


(5.1) 


where ^ ^ comparison of q with (3.12) shows that, as an operator, 

q = —id^ = —iX^. (5.2) 


The relation (4.11) implies the quantisation of the mechanical momentum according to 


.Vo 

p -)■ -i—+tAo^, 
or 


(5.3) 


where A is the magnetic monopole vector potential (4.6) 


^ If p were to vanish then (4.50|) forces q to vanish, and then (4.49) would imply E = 0, which is impossible. 
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Inserting (5.2) into (4.19) gives 


1 / 




V 


ip 


2V 


(5.4) 


which tnrns ont to be precisely the ganged Laplace operator (3.31). To see this, note that 


1^2 ^ id.yAi)(-idi + id^Ai) 


dr 




(6.6) 


where is the Laplace operator on Enclidean and we have nsed that, for the Dirac 

monopole ( |4.6 ), div^d = 0. In terms of spherical coordinates and (4.6) one checks that 

4 ) = + ^{dj + cot(3df^ + csc‘^(3dl), A~= ^ 


SO that 


df r 2 sin^ (3 


r 2 sin^ 13 ’ 

(6.6) 


IpI^ = —^dr{r^dr) —^(9^ + cot (3dp + csc^ (3d‘l^ — 2 cot /? esc (3d.yda + cot^ (3d^) 
= -^drir^dr) - + Xl), 

where we have nsed the relation, 

X^ + X 2 = d^ + cot I3di3 + cot^ (3d'^^ + csc^ (3d‘^ — 2 cot /? esc /3dad.y, 


(5.7) 


(6.8) 


which can be obtained from (3.12). Snbstitnting (5.7) into (5.4) shows that the qnantnm 
Hamiltonian Hp is the ganged Laplace operator (3.31), as claimed. 


Finally applying the qnantisation rnle to the angnlar momentnm J dehned in (4.22) we obtain 
the differential operator 


d ^ 

J = —if X — + i(f X A — f)dy 
dr 

( sin adi 3 + cot [3 cos a da — ^7 

— cos adp + cot j3 sin a da — I • 

-da 


\da 


(6.9) 


Transforming coordinates according to (3.10), one checks that, np a factor of i, the components 
are the vector helds Zi, Z 2 and Z 3 (3.14) generating the left-action of SU{2) on itself; 

J = iZ. (5.10) 


It follows that the sqnared total angnlar momentnm operator can be written in terms of the 
left- and right-generated vector helds on as 

.P = -(^2 ^ ^2 ^ ^ _^^2 ^ ^2 ^ ^ 2 ) ^ 
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5.2 Separating variables 


For fixed r, the angular part of the quantum Hamiltonian (3.31) is akin to the Hamiltonian of a 
symmetric rigid body coupled to a gauge held. In that context, the operators iZj are interpreted 
as ‘space-hxed’ angular momentum components and the operators iXj as ‘body-hxed’ angular 
momentum components The quantum Hamiltonian Hp commutes with Zi, Z 2 , Z^ and 

with X 3 ; together, these generate the U{2) symmetry of TN space. 


To separate the radial from the angular dependence in the wavefunction, we therefore require a 
complete set of functions on SU{2) which diagonalise the commuting operators 
This is usually done in terms of Wigner functions of the Euler angles, but here we use the 
construction of the eigenfunctions as homogeneous polynomials of the complex coordinates Zi, Z 2 
and their complex conjugates given in [3]. As explained there, an irreducible representation 
of SU{2) can be given in terms of polynomials in zi, Z 21 Zi^ Z 2 that belong to the kernel of the 
differential operator □ = 4:{didi + c? 2 ^ 2 )- Combining this observation with the discussion in [16] 
we obtain a basis that satisfies this irreducibility condition: 


sm 


{j + m)\{j — m)\ 

(5.12) 


(i + ^)' _ (j-m)!(-l)-^ ^ j+rn-k k^j-s-k 

^ (j + m — k)\k\ (j — s — k)\{s — m + k)\ ^ ^ 1 2 > 


where ^ 

j G-Z+, s,m =-j,-j + l,...,j-l,j, (5.13) 

and k runs over the values so that the factorials are well defined. These functions are nor¬ 
malised and are clearly orthogonal since they are eigenfunctions of the Hermitian operators 
A 53 , (total angular momentum), iZ^ (angular momentum along the space-hxed 3-axis), iX^ 
(angular momentum along the body-hxed 3-axis) with eigenvalues 


=-i(i + tZ3YL^<nY>„ iX^L = sYi^. (6.14) 

They also satisfy 

A'+yi„ = -ilU -s)U + s + X.Y,l = -i[(j + s)(j - s + (5.15) 


where X± = Xi ±iX 2 , which shows that all the angular momentum eigenstates can be obtained 
from the holomorphic Yj^ or the anti-holomorphic by the repeated action of X_ or X+. 

We look for stationary states T of the form 


T(r, zi, Z 2 ) = R{r)Y^^^{zi, Z 2 ). 


Using 


+^2)ni 


in the stationary Schrodinger equation 


HU + 1 ) + AyL 




(5.16) 

(5.17) 

(5.18) 
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we obtain the radial equation 


~dr{r‘^dr) + 4i(i + 1) + 





R{r) = 0. 


(5.19) 


Before we study bound and scattering states in the following sections, we make two general 
observations. 


It follows from (5.2) that, when acting on the functions (5.12), the operator q has the eigenvalue 
—s. For later use, note that the classical bound (4.21) also holds in the quantum case, so that, 
in particular, for any eigenstate of Hp and q with eigenvalues E and —s, we have 


L^E > 



if ps > 0 
if ps < 0. 


(5.20) 


Finally, It is worth stressing that neither the space-hxed nor the body-hxed angular momen¬ 
tum operators discussed above are invariant under f/(l)-gauge transformations. The quantum 
numbers j, s and m are not gauge invariant either and therefore have to be interpreted with 
care. However, this is familiar in the context of the Schrodinger equation coupled to a magnetic 
held, particularly in the discussion of Landau levels for planar motion in a magnetic held. Even 
though the angular momentum operator is not gauge invariant in this context, the eigenvalues 
can be used to label degenerate energy eigenstates. This labelling is not gauge invariant, but 
physical quantities like the energy or the degeneracy of energy levels are. The role of the gauge 
choice in labelling degenerate states in Landau levels is discussed in detail in HZI, see also the 
book [T8] . 


5.3 Bound states 

The substitution of 


, s - 2\ 2 

_ '-'2 —i 2 


R{r) = Re ^u{r), k'^ = 


E, 


into the radial Schrodinger equation (5.19), reduces it to 

(Pu ,du r ^ r. 


where 

and 


^ = 2k'r, a = (j -I- 1 -I- A), b = 2j + 2, 


^ ~ ~Wl + PS- 2s^) 


(5.21) 

(5.22) 

(5.23) 

(5.24) 


The equation (5.22) is the conhuent hypergeometric equation [19]. The general solution which 
is regular at the origin is 

u = ARe~’'^'''M{a,b,z), (5.25) 
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where A is an arbitrary constant and M is Knmmer’s fnnction of the hrst kind. Sqnare inte- 
grability reqnires 

a =—u, i/ = 0,l,2,... and /c'G M"*". (5.26) 

Since j takes arbitrary half-integer positive valnes, the hrst condition is eqnivalent to 

n := —A = u + j + 1, 1 / = 0,1, 2,.... (5.27) 


In principle, n can take all half-integer valnes > 1, bnt the ranges of the qnantnm nnmbers j 
and n are related by 

n = j + l + v, z/= 0,1,2,.... (5.28) 


This reqnirement together with the expression (5.24) for A as well as L > 0 and fc' > 0 imply 

L^E>2s^ -ps. (5.29) 


On the other hand, the relation (5.21) between E and k' enforces 

^2 


f py 2 , P 

L E<[s--'j =s -ps+—. 


(5.30) 


There can only be bonnd states if these two ineqnalities can be simnltaneously satished, i.e., if 

^ (5.31) 


which is the qnantnm version of the condition (4.50) for bonnded orbits in the classical theory. 


Note that, if L were negative, the ineqnality (5.29) wonld have the opposite direction and there 
wonld be no condition on p. In that case we can set p = 0 and recover the bonnd states in the 
singnlar L = —2 TN space discussed in [13], which exist for any s ^ 0. As shown in [20], there 
are no bound states (and no bounded orbits) when L > 0 and p = 0. More generally, however, 
binding is always possible when p is sufficiently large. All these results confirm the qualitative 
discussion of the two-dimensional toy model in Sect. 


Solving (5.21) and (5.27) for E, we find 


-I 


2,2 PS 
-n + s —— 




2 , P 

^ + “T- 

4 


(5.32) 


Only the solution with the upper sign satisfies (5.29), and we write the resulting spectrum of 
bound state energies as 




2 PS , ^ . /^2 _ .2 , P ^ ^2 


s -—s2 -(-——n 

2 V 4 


n = |s| + 1, |s| 2, |s| + 3 


(5.33) 


The behaviour for large n is typical for Coulomb bound states 

2 


E 


L2 


q 2 _ 2l 

* 4 


4L2n2 




(5.34) 
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This formula for the bound state energy shows that, like in the toy model of Sect. 2, the bound 
state energies are relatively high when p and q have the same sign (so that the signs of p and s 
are opposite) but are lowered when the signs of p and q are opposite (and those of p and s the 
same). Note also that, in the limit p = 0 and for L = —2 our formula reduces to that obtained 
in [13] for the negative mass TN space. For a detailed comparison observe that in [13] only 
integer values of j and n were considered. 


The energy levels for fixed s and n have a large degeneracy, given by the sum over the dimension 
2j + 1 for allowed values of j. Recalling the constraint (5.28) and j > |s|, the degeneracy is 




i=bl 


2j + 1 = (2|s| + 1) + (2|s| + 3) + ... + (2n - 3) + (2n - 1) = 


(5.35) 


with indicating that we sum over integers if |s| is an integer and over half-odd integers if 
|s| is a half-odd integer. As we shall see in Sect. 6, the degeneracy can be understood in terms 
of a conserved Runge-Lenz vector. 


5.4 Scattering states 


Next we turn to solutions of the eigenvalue equation (5.18) which describe stationary scattering 
states. For the analysis of scattering it is convenient to use parabolic coordinates familiar from 
the treatment of Coulomb scattering. 


Assuming solutions of (5.18) of the form 

T = e-“^e*™“A(/l,r), 


(5.36) 


and recalling the formula (5.8) and we hnd that A has to satisfy the equation 

[dl + ‘^dr + \ cot pda] A -I- ^(—cot^ P — rrp csc^ P -f 2ms cot P esc P)A 

\ r a ^2 J ^2 




L2 V“ 2VJ 

Now introducing parabolic coordinates (^, p via 

^ = r(l-1 -cos/3), p = r(l — cos/3), 
and noting the inverse transformation 

^ + V n . n 

r = -, cos/3 = --, sm/3 = 


A + EVA = 0. (5.37) 


(5.38) 






(5.39) 


we think of A now as a function of ^ and p via the substitution (5.39) for r and /3. Then (5.37) 
becomes 

4 




+ pd^ + d^ + d^)A-^ + m^ + 2ms|^^ A 


■p- + f + EL 

Jj 1 j 


2A 


p 




(5.40) 
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Separating variables again via A = f{^)g{ri) we find 


4e 


2 s^ ps 


—9^/ + + s) + k ^ + 2 ^EL — + ^ ) — C — 0, 

—dig + -dr^g - -(m - s)^ + k‘^g + C = 0, 

g ' g g 

where C is a separation constant and 


e = E-^{s-^ 


L2 V“ 2 

These differential equations can be simplified further if we assume solutions of the form 

^ V |m —s| ikE , ^ ^ ^ Im + sl ikn 

/K) = (j(ri) = ri^e-^G{ri), 


(6.41) 


(5.42) 


(5.43) 


(5.44) 


and implement the change of variable 

zi = ikE, Z2 = ikp. (5.45) 

Doing so we deduce that both E and G satisfy the confluent hypergeometric equation 

(5.46) 


(PE ^dE 

Zl-TT + (oi - - ciiP = 0) 


dzl 

d^G 


dzi 

dG 


where 


\m — s 


Oi = 


^ 2 —rY 4 " (^2 ~ ^2)3- 0.2G — 0, 

(IZ 2 (lZ\ 


1 ic i f 2 s^ ps\ , 


|m + s| 1 ic 

CLo — - H“ ~ “1“ TV’I ^2 “ “h 5 “h 1. 

2 2 Ak ' ' 


We see from these relations that 


3 1 | I 1 | I -x 

Oi + 02 = 1 + -\m + s| + -\m — s| — zA, 


where 


A = 


2kT 


i^L^E + ps — 2 s^) . 


So, the scattering solution is of the form 


^ = e-“V™"e^h“e-^e-^M(ai, 61 , ikE)M{a 2 , 62 , ikg), 


(5.47) 

(5.48) 

(5.49) 

(5.50) 

(5.51) 

(5.52) 


where M is Kummer’s function of the first kind. This is formally the same as the two monopole 
scattering state found by Gibbons and Manton ra in the case p = 0 and L = —2. However, 
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in our case the constants k and A have an extra p-dependence and the length parameter L is 
positive. 

As in [13] one can compute the cross section by looking at the wave function with m = s and 

ai = 1, 

^ — z)''^''M{\s\ — i\,2\s\ + l,ik{r — z)), z = rcos(3. (5.53) 

Then the substitution of the asymptotic form of M(|s| — i\, 2|s| + l,ik{r — z)) for large \z\ as 
in na allows us to identify the scattered spherical wave, and to obtain the cross section 


da 

dVt 


s2 + A2 
4A;2 


csc"^ |. 


(5.54) 


Writing A in terms of k via (5.51), we hnally arrive at the cross section 

2" 


k da\ 


16 

16 


4s2 


+ 1 - 


+ 






k‘^L'^ Ak‘^L‘^) 

2 


CSC^ I 


2 \ 2 

. S \ P 

1 H-H- - — 

k^L^) Ak^L^ 


2s2 

2 - + 


p 


A 


PL2 4PL2/ 


CSC^ I. 


(5.55) 


The special cases s = 0 and p = 0 are interesting because the resulting cross sections 


da A 

^ ) 


(pP) 


“ 16 * 


and 


k da\ 


(0,s) 


{p/^f 


= Y 

16 V ^ W ) 


CSC^ I 


CSC^ I 


(5.56) 


(5.57) 


are mapped into each other under the exchange s -H- p/2 even though the general case (5.55) 
is not invariant under this exchange. In both special cases, the charge, energy and angular 
dependence approaches that of the Rutherford scattering cross section for electrically charged 
particles in the limit of large s (or p). 


6 Algebraic calculation of quantum bound states 

6.1 The Runge-Lenz operator and so(4) symmetry 

In 1926, Pauli computed the quantum spectrum of the Kepler problem by using the conservation 
of the Runge-Lenz vector m His method has since then been much explored and extended in 
various papers, see [22| for a reference which is particularly useful in the current context. More 
recently, it was used to compute bound states and scattering or the Laplace operator on TN 
space [15] and also for the Dirac operator on TN [TT] . We now use it to re-derive the spectrum 


of the gauged TN Hamiltonian (3.31) purely algebraically. 
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As always in quantising a theory, we need to be careful with ordering in the quantisation of 
classically conserved quantities. While there are no such ambiguities in the dehnition of the 
angular momentum operators, they do arise in dehning a quantum version of the Runge-Lenz 
vector. The quantum analogues of the canonical Poisson brackets (4.17), 


IPi.Pil = [Pj./Wl = -i9jf{r), 


(6.1) 


imply, for the quantum angular momentum operator (5.9), 




( 6 . 2 ) 


which is the quantum version of (4.23). This means that [Ji,Pj] 7 ^ 0, i ^ j and hence the 
order of J and p is important in the dehnition of the quantum version of the Runge-Lenz vector 
( [T^ . 

Noting that, classically, p x J = ^{p x J — J x p), one hnds that the quantum ordering 


M = l{p x j-jXp)-^ (L^H - 27 ) 


(6.3) 


ensures that the quantum commutation relations between J and M are. 




(6.4) 


in analogy to the classical Poisson brackets (4.24) and (4.26) respectively. Now we use (6.2) to 
rewrite the Runge-Lenz vector as 


M = p X J — ip 




(6.6) 


The second ambiguity has to do with the position of the factor f of the last term. The above 
choice guarantees that the quantum Runge-Lenz vector commutes with H. 

In order to obtain a Runge-Lenz vector which commutes with the gauged Hamiltonian Hp and 
still satishes the relations (6.4), it turns out that the addition of the term / (4.28), which worked 
in the classical case, also works in the quantum theory. The gauged quantum Runge-Lenz vector 
is therefore 


MP = pxJ-ip- — {L^Hp 
A lengthy calculation yields the commutators 

k . m'h [q+rf-H, 


2q^ - pq) . 


( 6 . 6 ) 


^ijkJkj 


(6.7) 
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which quantise the Poisson brackets (4.30). We also hnd the following operator identities: 


MP-J = 
MP-MP = 


J-MP 

Hp- 


JL 

2L 



{L^Hp 

1)1 


- - pq) , 

{J ■ J — 0 ’^ + !) + 


1 


{L^Hp-2q^-pq)\ 


( 6 . 8 ) 


Since the Hamiltonian Hp and the U{1) generator q commute with each other and all compo¬ 
nents of MP and J, we can £x their eigenvalues and study the commtutation relations of Mp 
and J in a fixed common eigenspace of Hp and q. Denoting the eigenvalues by, respectively, E 
and —s, and assuming the bound state energy range 

(6.9) 

we dehne the rescaled Runge-Lenz vector, 

MP = ^ ^=MP. ( 6 . 10 ) 

^jUs-W-E 

Together with the components of J, it satishes the so(4) commutation relations, 

[Ji,Jj\=ieijkJk, [Ji,MP] = ieijkMl, [MP, Mj] = ieijkJk- ( 6 . 11 ) 


6.2 Bound states revisited 


The bound state energies of Hp can now be derived from the isomorphism so(4) ~ su{2)(Bsu{2) 
and the standard representation theory of su{2). We introduce the commuting operators 

J^ = ^{J±MP), (6.12) 

and see that the two Casimirs 

4 = \i-J- J + MP-MP)± ^{MP ■ j + j ■ MP) (6.13) 


have eigenvalues j±{j± + 1), where j± are both non-negative half-integers. 
J = J+ -|- J_, it follows that the total angular momentum quantum number j 
lies in the range 


Moreover, since 
dehned in (5.14) 


\j+-J-\ < 3 < \3++3- 


(6.14) 


Since —j<s<j, we deduce that 


- \ j++j-\ <s< \j++j- 


(6.15) 
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In terms of the relations (6.8) read 


■ J = J -MP = 


s {L^E-2s^+ps) 

V 

2 ^ {L^E-2s^+psf 


'U^-W-E 

+ J ■ J = -I + 


4((s- |)2-L2E)' 


(6.16) 


Substituting these into ( |6.13| ) and replacing by the eigenvalues j±{j± + 1), we get two 
quadratic equations for the unknowrj^ 


n := 


E^E — 2s^ + ps 

2lJUs-W-E 


(6.17) 


namely 


+ 2sn + — 1 — 4j+(j+ + 1) = 0, 
— 2sn + — 1 — + 1) = 0. 

(6.18) 

The roots of the hrst equation are 


n = -s± (2j+ + 1), 

(6.19) 

and the roots for the second equation are 


n = s± (2j_ + 1). 

(6.20) 

Both equations have to be satished for some values of j+ and j_, but combining the upper sign 
in one with the lower sign in the other implies a value of s which is outside the range (6.15). 
Hence, there are only two possible solutions for n, one which is manifestly a half-integer > 1 

n = -s + (2j+ + 1) = s + (2j_ + 1) 

(6.21) 

and one which is manifestly a half-integer < — 1 


n = -s- (2j+ + 1) = s - (2j_ + 1). 

(6.22) 


Finally solving (6.17) for E we obtain again the solutions (5.32) previously obtained via square 


integrability arguments. However, we still have four possibilities in total: two choices of sign 


in (5.32) and two choices for n (positive or negative) and in this section we cannot assume the 


conditions (5.29) and (5.31) to resolve the ambiguity. In the four different cases the energy 


equation (5.32) takes the two possible forms 


E^E 


2 , ps 2 

" +Y = " 



El _ q2 

1 + ^^^- 1 




(6.23) 


“^Note that the definition of n here is consistent with (5.24) and (5.27) 
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We can eliminate the lower sign becanse it conflicts with the lower bonnd (5.20). To see this 
we re-write (6.23) with the lower sign as 


L^E = -ps - 2{n^ - s^) - 2n‘^\ 1 + 


21 _ c2 

J^ 

9 ' 


(6.24) 


showing that L'^E < —ps in this case. However, this is inconsistent with the energy ineqnality 


(5.20) and therefore rnled ont. 


Tnrning to the npper sign, we need to consider the two possible signs of n and check the 


consistency between (6.23) and (6.17). In the case n > 1, both (6.23) and (6.17) assign a 


positive sign to -|- y provided < p^/4. In that case we arrive at the previonsly 


derived energy spectrum (5.33) together with the condition (5.31) for bound states. However, 
n < —1 is also consistent provided > p^/4. We have not been able to eliminate this case 
using only the algebraic methods of this section. It seems that the consideration of the actual 


wavefunction (5.18) and integrability requirement (5.27) is needed to rule out n < —1 


Finally turning to the degeneracy of the energy levels, we see that the quantum numbers n and 
s are determined via 

n = ]++ ]- + l, s=j+-j-, (6.25) 

and that the degeneracy of the energy level with quantum numbers n and s is the dimension 
of the tensor product 0 j_ of the irreducible angular representations with spins and j_. 


(2j+ + l)(2j_ + l)=n^-s^ 


(6.26) 


reproducing and interpreting the degeneracy (5.35) of energy levels. 


In this section we have only studied the commutation relations of the angular momentum and 
Runge-Lenz vectors at energies satisfying E^E < (s —p/2)^ and corresponding to bound states. 
It is not difficult to modify our discussion for the case E^E > (s — p/2)^. The dynamical 
symmetry algebra of the angular momentum and of a suitable rescaled Runge-Lenz vector. 


analogous to (6.11), turns out to be isomorphic to the Lie algebra so(3) x of the Euclidean 
group when E^E = (s — p/2)^ and isomorphic to the Lie algebra so(3,1) of the Lorentz group 
when E^E > (s — p/2)^. In the following section we will see how these three cases can be 
understood from a unihed, geometrical point of view. 


7 Twistorial derivation of the gauged Runge-Lenz vector 

7.1 Twistors, SU{2,2) symmetry and moment maps 

It has been known for a while [23] that the usual Kepler problem can be regularised by em¬ 
bedding momentum three-space into the three-sphere by means of a stereographic projection. 
This gives a geometrical picture of the angular momentum and Runge-Lenz vectors as con¬ 
served quantities associated with symmetries of the round three-sphere. For a full geometrical 
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understanding of the dynamical symmetry of the Kepler problem, it is moreover convenient to 
think of it as the symplectic quotient of an eight dimensional phase space, and the dynami¬ 
cal symmetry algebra for the various energy regimes as subalgebras of so(4, 2), see [21] for a 
pedagogical review. 

It was shown in [25] that one can similarly interpret angular momentum and Runge-Lenz vectors 
of the (ungauged) TN motion as generators of a subalgebra of a su{2, 2) ~ so(4, 2) symmetry 
algebra acting on an eight-dimensional phase space of twistors. In this section we show how 
this story can be extended to the gauged TN dynamics. We begin with a brief review of the 
relevant notation. 

For our purposes, twistor space is T = (C^ x C^)\{0}, and a twistor 

z = (“) e T (7.1) 

is a pair of spinors u = and vr = . Twistor space T is endowed with a pairing 

[Z, Z) = TfiCUi -|- 7f2(U2 + CUiTTi -|- ('^•2) 


which can be written as the matrix product Z*Z where 

z- = ztA A=(° ;»). 


(7.3) 


is the conjugate spinor and we write again tq for the 2x2 identity matrix. The pairing is 
invariant under f/(2,2), but we are particularly interested in the Lie algebra of the subgroup 
SU{2,2). Following [21], we pick generators 'Jkl where the structure constants are purely 
imaginary, and again write r* for the Pauli matrices: 


7ofc 

7fe6 

756 


Tk 

0 

0 

-Tk 


0 

-Tk 

Tk 
0 


To 

0 


0 

-To 


7q 2^7^ 


Tk 0 

0 Tk 


705 = 


0 To 
-To 0 


706 — 
1 


Ikb — 


0 

2 I To 
0 Tk 
Tk 


To 

0 


{.Ilk — —^KL, K,L — 0, 


2 \ Tk 0 

,3,5,6; i,j,k = 1,2,3). 


(7.4) 


For us, two sub-Lie algebras will be important. The stabiliser Lie algebra of the generator yoe 
is the Lie algebra generated by Si := ^eijk'yjk and ti := 75 ^ and is isomorphic to so(4): 


[Sj, Sj] if^ijk^ki i^ijktki [^*1 'i^-ijkSk- 


(7.5) 


The stabiliser Lie algebra of the generator 705 is the Lie algebra generated by Sj and := 7*6 
and is isomorphic to so(3,1): 


[Sj, Sj] i^^ijkSki ['^*5 i^ijkTki i^ijkSk- 


(7.6) 


29 


(7.7) 


As discussed in |25n the space T has the U{2,2) invariant one-form 


whose exterior derivative 


e = lm{Z*dZ) = —XZ*JZ^ - Z^dZ*J, 

2i% 


n = d9 = -idZ* A dZ 


(7.8) 


is a symplectic form on T. The su{2, 2) generators jkl dehne vector helds on T whose moment 
maps are 

Jkl = Z*^klZ. (7.9) 

The diagonal t/(l) subgroup of f/(2,2) acts on T, preserving its symplectic structure. The 
moment map is ^Z*Z, and the symplectic quotient is the level set 


Tg = <z eT 


-Z*Z = q 
2 


quotiented by the diagonal U(l) action: 


M, = TJU(1). 


(7.10) 


(7.11) 


We now introduce coordinates on T which are particularly well adapted for describing this 
quotient. With the notation r = (ti,T 2 ,T 3 ), we parametrise the spinors vr and oj in terms of 
spherical coordinates (R,a,f3,'y) and P G M.^,q G M as 


vr 


= Vr 


(e 2 ("+')') cos I 


V 62*^" 7) I J ’ 


and 


cu = . r + —To ) TT. 


(7.12) 


(7.13) 


In order to compute the symplectic structure and moment maps in terms of R and P, we note 
that 

(7.14) 

where 


ttV = P, vr'^^rTT = P, 

P = (Xi, X 2 , X 3 ) = (Psin /3 cos a, Psin (3 sin a, Pcos/3), 


and 


Juj = Rp‘^ + ^, Jfuj = -2pxJ+(Rp‘^ + ^]R. 

R \ RJ 


q 


(7.15) 


(7.16) 


It is not difficult to check that, for hxed q, the twistor Z = satishes (7.10) and thus belongs 
to Tq. Moreover, the diagonal U(l) acts simply by shifting 7, so that the vectors P,R E 
which are independent of 7, are good coordinates on the quotient Aig. 


^Note that our sign conventions differ from those in [25] 
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The symplectic structure (7.8) induces a symplectic structure on A4g which can be expressed 


as 

n = dXiA dPi + ^eunXndXi A dXi. (7.17) 

The moment maps for 750 , yeo and the generators of their stabiliser Lie algebras can be written 
in terms of P, R as 

Z*-f,oZ = i(a;tcu - vrtTr) = ^ (i?p2 + | - p) , 

^*^06^ = = - i^RP^ P ’ 

Z*sZ = -(w^rTT + TT^ra;) = P x P + gP, 

+ P^ P, 

Z*rZ = - u^fu) = Px J-l (rP^ + ^-R)r. (7.18) 

2 2 \ R / 

We can summarise these formulae more neatly by introducing a variable k, which can take the 
values ±1. Then we write 

4 = i(pp2 + | + «:i?), (7.19) 

for the moment maps Z*'^^qZ and Z*'yQQZ with k = 1 and k = —1. We write 

J = P X P + gP, (7.20) 

for the vector or moment map Z*sZ, and hnally note that 

K = Px J- HpR (7.21) 

unihes the moment maps Z*tZ and Z*rZ for k = 1 and k = —1. 

It then follows from the general theory of moment maps (and can also be verihed directly) that 
the Poisson brackets of these moment maps are, up to factors of i, the commutators of the Lie 
algebra elements which enter the dehnition. In other words, the brackets are 

{Hp,R} = {Hp,K,} = 0 

{Jii Jj} ^ijkJki {Jii ^j\ ^ijkXki ^^^ijkJk- (7.22) 


Z*tz = --(vr^rTT + Oj^foj) = P x J - - fpp2 + 1 
2^ ^ 2 V P 


7.2 Mapping twistor space to the gauged Taub-NUT phase space 


Adapting the treatment of [2S], we shall show how the quotient M.q (7.11) can be mapped onto 
the symplectic quotient M.q (4.14) of the cotangent bundle of TN. The map between the phase 
spaces is not canonical, but it can be extended to a map on the evolution space, preserving the 
presymplectic (or Poincare-Cartan) two-form. 
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For a phase space A4 with symplectic structure u and Hamiltonian ff, the evolution space is 
X M, and the presymplectic two-form is a; -|- dll A dt, where t is a global (time) coordinate on 
M. The trajectories of the flow with Hamiltonian Id can be characterised as the vortex lines of 
u + dH Adt, i.e., the lines whose tangent lines are in the null space of u + dH Adt. Now consider 
extended phase spaces x M and Ad x M with symplectic structures and Hamiltonians (cu, H) 
on Ad and (a), H) on Ad, and time coordinates f on M and t on M. Then a map 

FiAdxM^AdxM (7.23) 

which satishes 

F*{uj + dH A di) = u + dH A dt (7.24) 

will map trajectories in the Hamiltonian system {AA,u,H) to trajectories in the Hamiltonian 
system For details and a pedagogical discussion of Hamiltonian trajectories as 

vortex lines of Poincare-Cartan structures see m- 

We should stress that, in contrast to the treatment of the ungauged case with negative L in 
|25j . and unlike in the usual Kepler problem, no regularisation is required in our case since our 
Hamiltonian is smooth and hnite on the entire phase space. For dehniteness we focus on the 
case 2 

L2i7p<(g + |), (7.25) 

which is relevant for bounded orbits. 


The two extended phase spaces we would like to map into each other are Ad^ x M with presym¬ 
plectic two-form 

q 

a = u + dHp A dt = dxi A dpi -\ - -eunXpdxi A dxi -|- dHp A dt, (7.26) 

and Adq X M with presymplectic two-form 

S = dXi A dPi + ^^eunXpdXi A dXi dHp A di. (7.27) 

The required map is most easily written down in terms of the coordinates (P, P, d) of Adg x M 
and the coordinates {p,f,t) on Adg x M. It takes the form 

P : Adq X M —)■ Adq X M, {f,p,t)^{R,P,i), (7.28) 


where 


R = 




-H 


VI 


P = 


t = 


p 






|r ■ p + 2 


1 




(7.29) 
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A lengthy calculation shows 


F*{dXi A dPi) = dxi A dpi + 


df+ r ■ dp) A dHp 


p* 




U^+W-H, 

^ dXj^ ^ —^’ii^Xyidxi A dxi^ 


and 


F* {dHp A dt) = dHp A dt — 


|(p ■ dr + r ■ dp) A dHp 


Combining these, we deduce 


F*E = a, 

as claimed. It follows that F maps solutions of the Hamilton equations 


dr 

dHp 

dp 

dHp 

dt 

dp ’ 

dt 

dr 

to solutions of the Hamilton equations 




dR 

dHp 

dP 

dHp 

di 

dP ’ 

di 

dR 


(7.30) 

(7.31) 

(7.32) 

(7.33) 

(7.34) 


Having seen that F maps trajectories to trajectories, albeit traversed at different rates, we 
conclude this section by showing how F relates observables. It is easy to check that F maps 
the angular momentum in Mg to the the angular momentum in Mg, i.e., the substitution of 


(7.29) into (7.20) gives the TN angular momentum 

J = r X p + qf. 


(7.35) 


The Hamiltonians and the Runge-Lenz generators, however, are related by pulling back with F 


together with rescaling. Substituting the expressions (7.29) into the Hamiltonian (7.19) with 
K = 1, one hnds the re-scaled gauged TN Hamiltonian, 


~ MHp-2q^-pq 

Mr, = — 




(7.36) 


and the substitution into (7.21) (again with k = 1) gives a rescaled Runge-Lenz vector (6.10), 

p X J - ^{HHp - 2q^ - pq)f 


K = 


ttM + W-h, 


(7.37) 
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8 Discussion and Conclusion 


Our results show that the inclusion of the magnetic held (3.16) in the discussion of dynamics 
on TN space is both mathematically natural and phenomenologically interesting. The gauge 
held preserves the geometrical U (2) symmetry of TN and here we showed that it also preserves 
the dynamical symmetry of the associated phase space. 

As suggested by our qualitative discussion of magnetic binding in Sect. 2, the form of the 
magnetic held ensures the existence of bounded orbits classically and of bound states quantum 
mechanically. As a result, TN dynamics with the magnetic held (3.16) combines aspects of 
two paradigmatic systems of classical and quantum mechanics - the Kepler problem and the 
Landau problem of a charged particle in a magnetic held - into one model that preserves and 
merges the most interesting features of both. 

In fact, one can understand many results of this paper qualitatively by thinking of gauged TN 
dynamics as a combination of a Landau problem in the cigar-shaped submanifolds of TN space, 
reviewed in Sect. 3.1 and a Kepler problem on the base. The magnetic held on the cigars 


acts as ‘magnetic plug’ which keeps trajectories in a bounded region and produces quantum 
bound states, provided the energy is sufficiently small and 1^1 < |p/2|. This last condition is 
precisely the condition for bounded trajectories in our toy model of Sect. 2. A further link 
with Landau states, explained in [T^, is that in the limit where the TN parameter e in (3.4) 
(which we set to 1 in our discussion) is taken to zero, the gauged TN problem actually becomes 
a four-dimensional Landau problem. 


Some of the details of our results deserve further comments. The spectrum of quantum bound 
states (5.33) depends both on the magnetic hux parameter p and the quantum number s. 
Importantly, it also depends on the relative sign of those two quantum numbers through the 
appearance of the product ps in the energy levels (5.33). This is a consequence of the breaking 
of the discrete symmetry (3.18). In view of the interpretation of s as a component of the body- 
hxed angular momentum, the dependence on energy levels on the alignment of this angular 
momentum with the magnetic flux may is essentially an (anomalous) Zeeman effect. 


In the interpretation of the TN manifold as a collective coordinate space of magnetic monopoles 
in mu^, the body-hxed angular momentum quantum number s is an electric charge. Adopt¬ 
ing that nomenclature, and recalling that the parameter p measure the magnetic flux of an 
external magnetic held, the equality of the purely magnetic scattering cross section and the 
purely electric cross section may be viewed as a manifestation of a magnetic-electric duality. 
As far as we are aware, this is not expected in this context, and it remains an open problem to 
hnd a simple explanation of it. 
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